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Sequential ballistic deposition (BD) with next-nearest-neighbor (NNN) interac- 
tions in a A^-column box is viewed a time-ordered product of x A^-matrices 
consisting of a single s/2^block which has a random position along the diago- 
nal. We relate the uniform BD growth with the diffusion in the symmetric space 
Hn = SL{N, R)/ SO{N). In particular, the distribution of the maximal height of 
a growing heap is connected with the distribution of the maximal distance for the 
diffusion process in H]\f. The coordinates of H]^ are interpreted as the coordinates 
of particles of the one-dimensional Toda chain. The group-theoretic structure of the 
system and links to some random matrix models are also discussed. 
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During recent years much effort has been devoted to theoretical analysis of properties of 
surfaces obtained by aggregation of particles. Several models describing various properties 
of clusters grown by different deposition processes deserved special attention. Among them 
we can mention the models of surfaces grown by Molecular Beam Epitaxy (MBE) (see, for 
example, Ref. Q), Polynuclear Growth (PNG) H-Q and by Ballistic Deposition (ED) 0- 



11| . All these models belong to the Kardar-Parisi-Zhang (KPZ) universality class [12|. Our 
attention is focused on ballistic deposition. In the simplest setting, one can assume that 
elementary units ("particles") follow ballistic trajectories with random initial positions and 
specified direction in space and adhere sequentially to a growing aggregate. Despite the ex- 
tremely transparent geometric formulation, the problem of stochastic growth in the ballistic 
deposition process remains one of the most challenging subjects in statistical mechanics. 

The definition of a standard one-dimensional ballistic deposition model with next- 
nearest-neighboring (NNN) interactions (for more details, see Refs. dMll) is as fol- 



lows. Consider a box divided in columns (of unit width each) indexed by an integer i 
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{i = 1,2,..., A^). The free conditions are assumed for left and right boundaries. At the initial 
time moment, {t = 0), the system is empty. Then, at each tick of the clock, t = 1,2, ...,T, 
we deposit an elementary cell ("particle") of unit height and width in a randomly chosen 
column, i. Suppose that the distribution of these nucleations on the set of columns is 
uniform. Define the height, hi{t), in the column i at time moment t. Assume now, that the 
cells in the nearest-neighboring columns interact in such a way that they can only touch 
each other by corners, but never by their vertical sides. This implies that after having 
deposited a particle to the column i, the height of this column is modified according to the 
following rule: 

h,{t + 1) = max[/i,_i(t), h,{t), hi+i{t)] + 1 . (1) 

If at the time moment t nothing is added to the column i, its height remains unchanged: 
hi{t + 1) = hi(t) . A set of deposited particles forms a heap (a pile) as it is depicted in Fig. 

mh. 

The available theoretical analysis of stochastic deposition focuses mainly on the properties 
of the enveloping surface hi{T). Recall the well-known scaling relation characterizing the 
surface of a growing aggregate: 



[VaTh{T)f^ - — 



iVl/2 



N 



.1/2 



N'"9 {r/N'") , (2) 



where the brackets denote an average over different realizations of the noise, the variable 
T = T/N is the averaged number of particles per one column, and the function g{u) is the 
rescaled variable u = t/N^^'^ has the following asymptotic behavior: g{u) ~ u^^^ for m <^ 1 
and g{u) ~ const for m ^ 1. Such a behavior is typical for many non- stationary growth 
problems. The exponent 1/3 has been found already in [l2| for the original Kardar-Parisi 



Zhang (KPZ) model and then observed in a variety of other growth models. In [16|, [17 
it was realized that, for flat initial conditions, the distribution of a rescaled PNG surface 
heig ht, T~^^^{h(i,T) — 2r), converges as r — )■ oo to the Tracy-Widom (TW) distribution 



18[, describing the edge states of random matrices belonging to the Gaussian Orthogonal 
Ensemble (GOE). Note that in the droplet geometry the statistics of growing surface is 
instead described by the edge states of the Gaussian Unitary Ensemble (GUE) [l7|. Height 
fluctuations in the KPZ universality class were measured in experiments, both in planar 



19H21| and in curved geometry in the electroconvection of nematic liquid crystals [20|, [21 
and a good quantitative agreement with TW distributions was found. 

The typical analytic analysis of a growth process described by the equation ([T]) consists 
of a few heuristic steps which are usually beyond a precise control. First of all, one sends 
the lattice spacing and time step to zero (with proper rescaling of all relevant variables), 
thus arriving at a continuous nonlinear partial differential equation of the type 

dthix,t) = iydlh{x,t) + g\dMx,t) \ + r]{x,t), (3) 

where ly and g are some constants and r]{x,t) is a Gaussian white noise, which replaces the 
telegraph noise in the initial equation ([1]). The next step is to replace the nonanalytic non- 
linearity \dxh(x, t) \ by the analytic one {dxh{x, t))^ This replacement is based on numerical 
evidence |22| and solid mathematical estimates j23| which show that the universality class 
of KPZ equation is insensitive to the value of a (1 < a < 2) for the nonlinearity of the type 
\dxh{x,t)\". 



4 A. Gorsky, S.K. Nechaev, R. Santachiara, G. Schehr 



Another point of view on ballistic deposition process has been developed recently in 24j . 
where it has been demonstrated that the discrete stochastic equation describing the BD 
process in ([T]) can be naturally represented in terms of a "dynamic programming" language 
associated with the so-called Bellman equation 25|]. This "dynamic programming" point of 
view allowed for a systematic translation of evolution of growing BD clusters in the NNN BD 
model into the language of maximizers and shocks in discrete equations of the Burgers or 
Hamilton-Jacobi type. This approach allowed to determine the exponents characterizing: i) 
the decay of the number density of clusters, ii) the wandering of the inter-cluster interface, 
and iii) the tail of the cluster mass distribution. 

In the present paper we consider the ballistic deposition process from a completely dif- 
ferent point of view, linking it to the diffusion in some symmetric space. The goal of our 
work is two-fold: 

1. On one hand, we conjecture that the uniform NNN ballistic growth in a scaling limit 
can be interpreted as multiplication of random matrices of type ( [TT|) . which, in turn is linked 
to the diffusion in the symmetric space Hjy = SL{N, R)/ SO{N). In particular, the distri- 
bution of the maximal height of a growing heap is connected with the distribution of the 
maximal distance for the diffusion process in Hj\f; 

2. On the other hand, we demonstrate that the coordinates of if at can be interpreted as 
the coordinates of particles of the one-dimensional Toda chain, which itself can be regarded 
as Inozemtsev limit of Calogero-type systems describing the transmission in random wires. 

The paper is structured as follows. In the Section [IT] we remind the concept of "local 
groups" and demonstrate that the BD process can be described in terms of matrix multipli- 
cations; in the Section IIIII we describe the NNN BD process as diffusion in an appropriate 
symmetric space, derive the connection with the one-dimensional Toda system and show 
that BD can be described by the Inosemtsev limit of hyperbolic Calogero-type systems; in 
the Section [IV] we interpret some group-theoretic facts concerning Toda chain in terms of 
ballistic deposition; the Conclusion summarizes our results. The two Appendices outline: a) 
the connection between heights and eigenvalues, and b) the link between ID Toda system 
and rational solutions of Painleve II. 



II. BALLISTIC DEPOSITION AND PRODUCTS OF RANDOM MATRICES 



A. Ballistic growth and directed polymers in random environment 



Consider the set of matrices {gi, ...,gN}, where the matrix Qi has the following form: 

/I \ 



1 
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u 
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row z. 



(4) 



1/ 
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The deposition event in the column i means the apphcation of the generator gi. Rising a 
heap by sequential addition of particles with NNN interactions is translated into a sequential 
multiplication of matrices of type Thus, the whole heap is represented by the time- 
ordered product 



V{T,u) 



T 
i=l 



(5) 



Each element of the matrix V{T, u) is a polynomial of the variable u. Take a vector a(t = 
0) = (ai, ...,aAr) where ai{i = 1, ...,N) are strictly positive values. The set of local heights 
h(T) = {hi{T), /i7v(T)) at time T after deposition event can be extracted as follows: 



h(T) 



lim 



\n[V{T,u) a{t = 0)] 



(6) 



The definition of heights in is consistent with the updating rules ([T]) for the NNN- 

ballistic deposition process. Namely, write the recursion for the updating dynamics of the 
vector a(t): 

\ / ai{t) \ 



/ai{t + l) \ 


/I 




aj^i{t + 1) 
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a,{t + l) 




u u u 


aj+iit + 1) 
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\a7v(t+l) / 


v 





aj_i(t) 
a jit) 
aj+i{t) 



(7) 



V \aN{t) I 

\iit = i (i.e. we drop a particle in the column i ar time t), then 

ttiit + 1) = u ai_i{t) + u ttiit) + u ai^iit). 

Otherwise aj(t + 1) = aj(t). Substituting the ansatz ai{t) = e^''*'^*^ into 
deposition event in the column i: 

Kit + 1) = ^In [e^'^-^W + e'^^'W + e'^^'+^W] + 1. 



(8) 

we get for 
(9) 



For all other columns j ^ i the height remains unchanged: hj{t + 1) = hj{t) [j = 1, A^). 
The initial configuration by definition is h(T) = (0, 0). Taking in Q the limit /3 — oo, we 
recover ([1]). Note, that iQ resembles the ultra-discretization of stochastic growth equations 



The quantity aj(t) can be interpreted as a finite temperature partition function of a 
directed polymer in a random environment, /3 = logu being essentially the inverse temper- 
ature. To demonstrate this consider the random lattice depicted in the Fig. [H The rules 
of the deposition process mean that at each "time" t a single site (an open dot) is chosen: 
an energy — e is assigned to this site and it is connected to three neighboring ancestors at 
time t — 1. To the other sites (black dots) an energy is assigned and they are connected 
(vertically) to a single site. On this random lattice we then consider a configuration of a 
directed polymer, shown by a bold line. 
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t=0 



t= 1 



t = 2 



t=3 



t = 4 




E(9) = 0; E{O) = -E<0 

Figure 1. Directed polymer model in a disordered environment. Open dots designate the randomly 
placed sites with energy — e. 



The energy assigned to a polymer configuration is the sum of site energies through which 



the directed polymer is passing. If we denote by e 



u. 



then the quantity ai{t) evolving 



according to ([7]) gives exactly the partition function of a polymer ending in site i at time 
t. By this correspondence a precise connection between the BD a directed polymer in a 
random medium is established (see, for instance, 27|] and references therein). 



B. Ballistic deposition and random walk on local groups 



We have seen that the ballistic deposition process has a natural description in terms of 
multiplication of random matrices with special "local" structure. The established connection 
is a particular example of a uniform Markov process on a "local group" j28l |. 

Definition 1. Local group (semi-group) Lj^ (L%) with N generators gi,...,gjyi is a group 
(semi-group) with the commutation relations 

9k9m = 9m9k, V|/c-m|>2, (10) 

and might have additional relations R inside the elementary blocks 9k,9k±i {k = 1, ...N). If 
the relations R are the same for all k, then L is called a local stationary group l2al . 

Plenty of important groups, semi-groups and algebras are of type of local groups: the 
Coxeter groups, Hecke algebra, braid groups, Temperley-Lieb algebra, etc. A locally free 
group F with M generators, being in some sense the "simplest" nontrivial example of a 
local group, can be considered as a universal object in the manifold of all local groups with 
a given number of generators. 

Definition 2. Locally free group (semi-group) (F^) with N generators gi,...,gN is a 
group (semi-group), determined only by the relation / f70|) . Each pair of neighboring genera- 
tors 9k,9k±i produces a free subgroup (sub-semigroup) of Fn (F^) \2!A }. 

Let us define the numbered heap as a heap whose cells are enumerated by sequential "ticks 
of time", t. Deposition of a particle in a column i means the application of a generator gi. 
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There is a bijection between numbered heaps of T sequentially deposited cells and words of 
length T in such that a given numbered heap is uniquely represented by a word in 
(and vice versa) j29|. The key point of this bijection is the special way of representing words 
in the locally free semigroup (group). Namely, all words are written in the so-called "normal 
order form" implying that the generators with smaller indices are pushed as left as possible 
when it is allowed by commutation relations ffTOj) . For example, the word associated with 
subsequent deposition of all 13 particles in Fig. [2] is: W = gs ga gi g^ gi g2 g5 93 Qi 95 93 96 92- 
The same word W written in the normal order form W = gigig3g2 9194,93 93 92 96 95 95 96, 
is in bijection with the numbered heap (Fig. ^) obtained by "contraction" of the sequential 
deposition (Fig. [2^). The statistical properties of locally free semigroups and groups are 



investigated in details in [29|. Thus, the uniform growth of a random heap on a set of A'" 
columns can be considered as a uniform Markov chain on a semigroup F^ with generators. 
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(b) 



1 2 3 4 5 6 



Figure 2. (a) The sequential ballistic deposition; (b) the heap obtained by contraction procedure 
using commutation relations of {N = 6) . 



It should be noted that the concept, equivalent to the locally free semi-group F^, ap- 
peared earlier sO] in the investigations of combinatorial properties of substitutions of se- 
quences and so called "partially commutative monoids" j3l|- Especially productive is the 
geometrical interpretation of monoids in the form of a "heap", offered by X.G. Viennot 
30| . It is worth mentioning that the concept of a monoid-heap appears in connection with 



various aspects of directed growth, parallel computing and complexity algorithms |32|. 

Definition 3. Take a locally free semi-group F^ with fixed finite number of generators 
{gi, gi\i} . Consider the uniform distribution on the set {gi, gjsf} . Construct the (left- 
handed) Markov chain on F^: a(t = 0) = {oi, oat}, a(t) = gi^a.{t — 1), where it = 
{1,...,A^} is uniformly distributed on the set of indices {!,... ,A^}. 

Analogously with locally free semi-group, we can define the Markov chain on a locally free 
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group. The geometric interpretation of random walk on locally free group is very straight- 
forward. Consider two types of particles: black and white, and identify the generators Qi and 



(yfj respectively with deposition of black and white particle in the column i 



TV). 



The elementary depositions are depicted in Fig. [3|, where (a) corresponds to the deposition 
9x9292 ^ = 9i (since 9i9~^ = 1), and the configuration shown in (b) cannot be reduced. 



(a) 



(b) 



Figure 3. Deposition rules for colored heap (annihilation is allowed). 



C. "Soft deposition" 



Deposition process viewed as a random walk on a locally free semigroup (or group) can 
be considered from the point of view of diffusion in symmetric spaces. Let us demonstrate 
this on two simple examples: sequential deposition in two and three columns. 

Deposition in two columns corresponds to a Markov chain on a locally free semigroup 
with two generators, 91,92 can be represented by a Cayley tree shown in Fig. |l]a, since 



for N = 
example 



2 one has F2 



F2 , where Fg is a free semigroup with two generators (see, for 



33| for definitions). Deposition in three columns can be treated similarly as a 
Markov chain on the Cayley graph associated with locally free semigroup F^: {91,92,93} 
and the commutation relation gig^ = g^gi - see Fig. |4]b. Certainly, the same correspondence 
can be established between the ballistic deposition with annihilation (as shown in Fig. [3]) 
and random walk on the full locally free group F/y. 




(a) 




Figure 4. Ballistic deposition in two (a) and three (b) columns labelled by states of Cayley graphs 
of corresponding locally free semigroups. 



We shall also consider a random walk on the full group Fj^, i.e. the uniform growth in 
a two-color heap in which direct deposition of white particle on top of black one (and vice 
versa) leads to their mutual annihilation. 
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The key problem connecting the A^-column deposition and diffusion in symmetric space 
is the following. For the multi-column deposition we need to find a metric space, H]sr in 
which the locally free group F/v with generators acts as a group of isometrics, 
which the Cayley graph Cn of Fn can be isometrically embedded. 



I.e. m 



It is well known that any regular Cayley tree, C, as an exponentially growing structure, 
cannot be isometrically embedded in an Euclidean plane. The embedding of a Cayley tree 
C into the metric space is called "isometric" if C covers that space, preserving all angles 
and distances. For example, the rectangular lattice isometrically covers the Euclidean plane 
E = {x, y} with the fiat metric (is^ = dx^ + dy"^. In the same way, the Cayley tree C 
isometrically covers the surface of constant negative curvature (the Lobachevsky plane) H. 
One of possible representations of H is the Poincare model, the upper half-plane Im z > 



of the complex plane z = x + iy endowed with the metric ds'^ 



of constant negative 



curvature. Thus, the locally free group F2 = T2 is the group of isometrics of Lobachevsky 
plane. Since the Lobachevsky plane is the coset space H2 = SL{2, R) / 5*0(2), we can consider 
the two-column ballistic deposition (with possible cancellation) as a random walk in H2. 

Similarly, the locally free group F3 generated by gi,g2,g3 with the commutation rela- 
tion QiQ^ = g^Qi acts as a group of isometrics of the Lobachevsky space with the metric 
ds'^ = '^^ ^'^y2 ^'^^ • Thus, the deposition in three columns has an important difference with 
respect to the two-columnar growth. Namely, the steps along the coordinates x and z in 
Lobachevsky space commute which is straightforwardly translated into the commutation 
relation (71(73 = (73(71 in the group F3. Hence, the three-columnar growth corresponds to the 
symmetric random walk in Lobachevsky space = SL{3, R)/ S0{3). The consideration of 
the corresponding diffusion equations in curved geometry is the subject of the next Section. 

To exploit the power of algebraic methods, we suggest the following smoothed description 
of A^-columnar ballistic growth. 

Definition 4. We define the "soft" deposition generator, gi, as follows: 



/I 



\ 



u 



1/ 



u. 



it) 




At) At) _ At) (t) 



where the block is an SL{2, R) matrix. For definiteness, suppose that the entries 
af\bf\cf\df^ are random variables (up to the condition det uj;*^ = ^) distributed in a 
compact support with some measure P{u). 

Note that this "smoothed" description captures all principal properties of the growing 
heap which we would like to analyze: 

i) gi{'ti)gj(t2) = gj(ti)gi{t2) for \i—j\ > 2 and t2 > ti, meaning that sequential deposition 
in distant (i.e. non-nearest-neighboring) columns commutes; 

ii) gi(ti)gj{t2) 7^ gj{ti)gi{t2) for \i—i\ = 1 and ^2 > ^i, meaning that sequential deposition 
in nearest-neighboring columns does not commute. 
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The group G : {gi, ^at} falls under the Definition 1 of local groups. The advantage of 
the redefinition of deposition generators Qi — cji [compare (jlj) and ffTTj) ] will appear below. 
Equipped with these definitions, we formulate the following conjecture: 

Conjecture. The Cayley graphs associated to the locally free group generated by gi 's on 
one hand, and to the local group generated by g^ 's on the other hand, can both be embedded 
isometrically in the same metric space H^- 

The metric space if at is the symmetric space if at = SL{N, R)/ SO{N). One can then 
use the Lie group structure of this symmetric space to define the Beltrami-Laplace operator 
which coincides with the Casimir operator C2. This would allow us to describe the symmetric 
Markov chain on the locally free group Fn with generators by the uniform diffusion process 
in the symmetric space Hn, where the group Fn acts as a group of isometrics. 

The connection between the columnar heights, hj{t) (defined in ([H])), and current coordi- 
nates fij{t) describing the diffusion process in the symmetric space Hn, is as follows. Denote 
by Vt the time-ordered product of t generators g in (ITT]) . Let {Ai(t), AAr(t)} be the set of 
eigenvalues of the matrix Vf, define fij(t) = In Xj(t) (j = 1, ...,N). Then, one has for t ^ 1: 

hraa.^{t) = 'J^ Atmax(^) (12) 

where /imax(^) = max{/ii(t), /lAr(t)}, /imax(^) = max{/ii(t), ...,/iAr(t)}, and 'Jn is some 
constant which depends only on the matrix size, N. Our numerical computations show that 
7Ar converges rapidly to its limiting value 7 ~ 0.8 as increases - see Appendix |Al We also 
describe in Appendix |X] the relation (for N = 2) between the number of steps, n, from the 
origin along a Cayley tree embedded into ii2 = SL{2, R)/ S0{2), and a hyperbolic distance 
d in H2. 



III. DIFFUSION IN SYMMETRIC SPACES 

A. Sequential deposition in two columns and a diffusion in if = SL{2, R)/ SO{2) 

For N = 2 the sequential deposition process is viewed as a multiplication of SL{2, R)- 
random matrices. The equation for the stochastic evolution of the vector r(t) = {x{t),y(t)) 
is 

x(t + l)\ /aW 6W \ /x{t) 
y{t + l)) ^ i cW rfW j [y{t) 



(13) 



which can be conveniently rewritten, introducing the Ricatti variable p{t) = x{t)/y{t), as: 

ol'l + 6<')p(«) 



The stochastic equation ( IT^ was the subject of the pioneering work [SJ] where it has been ex- 
plicitly shown that the limiting distribution function Wt{p) for T ^ 1 and af\bf\cf\df'' 
randomly distributed on a compact support with one and the same distribution for all 
t = 1, ...,T, satisfies the diffusion equation in the Lobachevsky plane (the two-dimensional 
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Riemann surface of constant negative curvature). Recall that the Lobachevsky plane coin- 
cides with the symmetric space H2 = SL{2, R)/ S0{2). So, we could immediately write the 
heat equation in H2 for the distribution function P(x, t), where t is the "time" and x are 
the coordinates of a Riemann manifold: 



dtP{^, t) = DAblP(x, t)- P{^, t = 0)= 5(x) (15) 

where 



9 / d 



is the Beltrami-Laplace operator and Qik is the metric tensor of the manifold; g = det Qik- 
The function P(x, t) satisfies the normalization condition J y/gW{'K,t)d:K = 1. 

It seems to be instructive to outline the derivation of (IT^ using the standard technique 
of extracting the Lie algebra generators from the appropriate matrix decomposition and 
writing the corresponding Casimir operator. Let us note that the proof of convergence and 
existence of limiting distribution associated with the stochastic equation f|T^ was a subject 
of the work [35| . 



The SI2 Lie algebra is defined by the relations [36 

[Xi, X2] = X3; [Xi, X3] = -2X1; [X2, X3] = 2X2, (17) 
where the generators Xi,X2,X3 have the 2x2 matrix representation 

X.-{ll): X.^(-); .3^(;_°,). (18) 
As we shall see later, the geometry of the system presumes the use of Gauss matrix 



decomposition [36|. So, we represent the S'L(2, i?)-matrix U as: 



In the differential form (in the right representation) the s/2-generators Xi,X2,X3 read: 

X, = d,; X2 = e'^^^dy±xd^-x'd^; = -2xd, ± d^. (20) 

The Casimir C2, which coincides with the Beltrami-Laplace operator Abl, can be expressed 
as follows 

C2 = Abl = X1X2 + X2X1 + ^X| = ±d^ + 2e^''^d,dy. (21) 
Thus, the diffusion equation for the distribution function Vr(/i, x, j/, t) has the following form 

^tW{^i, X, y,t) = D ±d, + 2e^^^d^di^ W{y., x, t), (22) 

where D is the diffusion coefficient. Let us make a comment concerning two options to 
choose the sign in equations f|T9l) - fl22|) . It is obvious that in f fT9|) the sign can be chosen 
arbitrary, however in fl2(I]) -f l2^ the sign should be consistent with correct physical behavior 
of the system at large times. The simplest way to chose the correct sign is to compare the 
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asymptotic behavior of W with the hmiting distribution of large number of noncommut- 
ing matrices. The classical paper by H. Fiirstenberg 37| states, that under some natural 
conditions typical products of T noncommuting matrices [SL{2, R) belongs to the class un- 
der consideration] increase exponentially with T. This selects C2 = — + 2e^^dxdy. 
Seeking now the solution of f l22|) in the form x, y, t) = e"*+''^Z(/i, x, t) and choosing 
a - 



-y, 6 = 1 we arrive at the Liouville-type equation for the function U{fj., x,y,t): 



D 



dtZifi, X, y,t) = - {dl + Ae'^d^dy) x, y, t). 



(23) 



In what follows we are interested in the //-dependence ("radial part") of the function U 
only. In this case the "angular" dependence on x and y can be eliminated by the reduction 
conditions 3^, which imply the replacement dxdy — )■ = const. Thus we arrive at the 
equation 

d,Z{f,, t) = ^ {dl + Ax'e'^^) Z{fi, t), (24) 

subject to the initial condition Z(/i,0) = S{fi). The solution of (124|) is well known and can 
be written explicitly, for example, as in jsoj 



Z(/i,t) 



dktlJk{fJ'Wk{0)e 



= 2^/-^ smh -=i^2,,/^ ( 4y — e'^ 



(25) 



where Ka{z) is the Macdonald function. For generic forms of solution of (12^ see [40 . 



B. Sequential deposition in three columns and diffusion in H = SL{3, R)/ S0{3) 



As we have seen in Introduction, for = 3 the sequential deposition process is not 
reduced to a random walk in H2 and the situation is more involved. If one considers two 
sequential deposition events, we can encounter with equal probability all configurations 
schematically depicted below. Namely, besides the former SL(2, R) matrix multiplication 





sU 



(26) 



we have also the matrix multiplication of two "shifted" copies SL{2, R)+^ SL{2, R)^, where 
the notation +(— ) designates upper (lower) s/2-block 




(27) 



It can be easily demonstrated that the multiplications of type fl2B|l - fl271) do generate the 
whole group SL{3,R). We use the Lie algebra representation of two "shifted" algebras SI2 
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and check that all s/3 generators can be obtained only from the multiphcations of type 
f l26p -( P7|l . Using s/2-representation f llSp we construct the set of generators of two "shifted" 
copies of SI2: 



1 0> 




/ > 




/ 1 





j ; ^2 = 
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> -1 



(28) 



These generators satisfy the commutation relations of s/3 algebra: 

[Xi, X2] = X3; [Xi, Xs] = -2X1; [X2, X3] = 2X2; 

[Fi, 1^2] =1^3; [11,13] = -2ri; [r2,l3] = 2F2; . (29) 

[Xi, Fa] = Xi; [X2, F3] = -^2; [Fi, X3] = Fi; [F2, X,] = -Y^ 

However the list ( 12 8 p does not enumerates all generators of s/3: two generators of s/3, X4 
and F4, where 












^0 







: 



































X4 = 000 ; F4 = 000|. (30) 

yoooy \ioo 

are still missing. They can be straightforwardly recovered as 

X4 = [Xi, Fi] = XiFi; F4 = [1^2, X2] = ^2^2 . (31) 

Thus, by using the set of two triplets ((Xi, X2, X3), (Fi, F2, 1^3)) (and their commutation 
relations) we generate the whole group SL{3,R). In other words, the sequential matrix 
multiplication of type ( 12B]) - (127|) generates the group SL{3,R). Qualitatively, this can be 
seen by writing the square of the metric element of two copies of Lobachevsky planes, dsf 2, 
in the upper half-plane representation in subspaces (x, y) and {y, z) of the three-dimensional 
half-space {x,y > 0,z) (Siegel space): 

dsj = ^^2^^ for H+ = SL{2,R)+/SO{2) + 

y (32) 

dsl = ^^^^ for = SL{2,R)^/SO{2)_ 

The square of the metric element, ds"^ = ds"^ + ds'^ in the Siegel space {x,y > 0,z) thus 
reads 

ds^ = d-' + ^df + d^\ (33) 

By rescaling y — )■ y/ we arrive at the standard Riemann metric of the Lobachevsky space 
H3 = SL{3,R)/SO{3), 

d~s^ = '-ds^ = d-' + dy^ + dz^_ ^3^^ 



Equation fl51|) allows us to see that the subspaces {x, y) and (y, z) are two copies of 
Lobachevsky planes: and H2 ■ 
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One expects that after a sufficient number of depositions, the Haar measure is smeared 
over the whole group SL{3, R) such that we arrive at the random walk in the symmetric 
space H-s = SL{3, R)/ S0{3). The diffusion equation associated with (12^ - (I27|) is governed 
by the Beltrami-Laplace operator Abl- The most straightforward way of constructing it 
consists in writing the Casimir C2 for infinitesimal transformations on SL{3, R)/S0{3) (see, 
for example, [38|). Again using the Gauss decomposition, we can decompose the generic 
matrix U G SL{3, R) as follows: 





/ ail 012 oi3 \ 




/ 1 
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/ 1 Xl X2 \ 


u = 


0-21 CL22 023 




yi 
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1 X3 




\ 031 0,32 033 / 
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1/ 
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e~^" / 


\0 1 / 



(35) 



The infinitesimal generators of the group SL{3, R) in the right representation form are as 
follows 



X2 
X3 

Yi 
Y3 



Xi9^i - x\d^^ - XiX2d^2 - (X2 - XiX3)(9, 



9x2 

d^.^ + xid^^ 

xid.^, - 2x3(9^3 - X2d^2 + (9^2 

p/X2-2/xi 



Ml 



^3^X3 



(36) 



X3dy, - e'^^-'^^'-'xidy^ + + ^^3,^gM2-2Mi) ^ 3,^^^ 



Ml 



(X1X3 - X2)<9^2 - xiX2<9xi + a;3(a;iX3 - X2)d^., - xld^^ 



The Casimir C2 on SL{3^ R) has the following form: 

C2 = X1X2 + X2X1 + ^1^2 + Y2Y1 + + ^4X4 + I (X| + Yi 



X3Y3) . (37) 



Writing C2 explicitly in terms of the generators given in Eq. fl36|) . we get: 

IC2 = IAbl = I {dl + dl^ + d,,d,,) + d,, + d. 



Ml 

p-Mi-M2^ n 



■'M2 



-2^1+^2 



{dx^dy^ + X3d.j,^dy^ + y3dx^dy^ + X3y3d.j,^d, 



y2, 



(38) 

Let us introduce the "radial" function Z{jjLi, jjL2,t) such that iy(/ii, /i2, t) = e ^* (M1+M2) 
xZ(/ii, /i2, t) and apply to fl38l) the reduction scheme |38| as in fl2^ fixing = dy^ = 0. 
Denoting d^^dy^ = k\ and dx^dy^ = x^, where the constants xi and X2 do not depend on 
the radial coordinates /ii and /i2, we get the diffusion equation for the radial part Z: 



D 



9,Z(/ii, /i2, t) = - [dl + + d,,d,, + 3^^^--"^- + 3x2 



-2/.ti+/12 



] Z(/ii,/i2,t). (39) 



The symmetry of the problem allows us to set in the equation fl39l) xi 
corresponding invariant Riemann metric ds^ is 



X2 



X. The 



92 

Ml 



5mi^m2 + dl + e'^^^-^^-dxA. + e-^^-+^^-dx,dy, + X3y3e^^+^^dxA 



(40) 



^yi 
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As we shall see in the next Section, the diffusion equation f l39p coincides with the diffusion 
equation constructed with the one-dimensional Toda Hamiltonian for two particles in the 
Chevalley basis 4l| with eliminated motion of the center of mass of the system. Note that 
the fact that we obtain the Casimir C2 in f l38|) in the center-of-mass frame is the consequence 
of working with SL[2,R) matrices with determinant 1. 



C. Sequential deposition in columns and one dimensional Toda chain 



The generic case > 3 can be constructed by considering the algebra of shifted copies 
of local Lie algebras s/2. The Markov chain on the uniform distribution of A^ elementary 5/2^ 
blocks results in growth of the random band-like matrix schematically depicted as follows: 



\ 



■I 



/I 



\ 



■J 



/I 



sU 



( 



t matrices 



(41) 

In order to capture the physics of the ballistic growth process it is instructive to begin 
with non-rigorous, essentially intuitive, construction. Suppose for a moment that A^ is odd 
and the deposition events are such that the positions of s/2^blocks either coincide, or do not 
overlap at all. One particular example is shown below: 



\ 



\ 



■I 



/I 



v 



\ 



J 



/I 



V 



\ 



J 



(42) 



In this case we can easily associate the growth process with the radial diffusion on N jl 
noninteracting copies of Lobachevsky plane H2 = SL{2, R)/ S0{2). By fl24p we can easily 
construct the radial part of the Hamiltonian of A^ noninteracting particles 



N 



(43) 



where pj = is the momentum of the particle j (recall that in the deposition process j is 



the column number), while /ij is associated with the height in the column j. 

The simplest way to incorporate the interactions between nearest-neighboring columns 
which flatten the growing surface, making it less rough, is to replace the Hamiltonian H^^^ 
of noninteracting particles by the following one 



TV 



H 



Pi 



KC 



') ; fiN+i = 0. 



(44) 
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We obtain this Hamiltonian by direct computations for = 2, 3 and bring also in Section 
HV] group theoretic arguments in favor of fl44p for arbitrary A^. 

First of all, let us demonstrate that the diffusion equation ( 139|) obtained by Gauss de- 
composition of the group SL{3, R) do coincide with the two-particle Toda Hamiltonian 



H{2) =pj+pl + Kef"'- 



'M2 



under the special choice of the coordinates. Actually, in the frame of the center-of-mass for 
the A^-body system, the variables fij [j = 1, ...,N) are transformed as follows (il ] 



yUl /il + 
1 



■ l^j + /^i + TTTT 5 2 < J < A/; fiN+i -l^N + 



N +1 

Pn+1 



N + 1 

N j-l 

i=j 



i=l 



N+1 
+ PN+l A<J<N 



N + 1 



N 



1 ^ 

1=1 



ipi + Pn+1 ■ 



(45) 

The explicit form of the A^-body Toda Hamiltonian in the coordinates ( 145|) is given in 
42|. For our purposes it is sufficient to check the A^ = 2-case only. The corresponding result 
is: 

H{2) ^\{pI+pI+ P1P2) + [e^^"''^^ + e-2'^i+'^2) , (46) 



which after the substitution pj 



^ gives (for x\ = h\ = k) the diffusion equation fl39l) . 



D. Toda chain and scaling limit of hyperbolic Calogero— type systems 

The Markov chain corresponding to the uniform sequential deposition of elementary s/2- 
blocks results in a growth of a random band-like matrix Vj^ii) [as shown schematically in 
f l4T|) ] with a time-dependent bandwidth, h{t\ and Gaussian distribution of matrix elements 
inside the band. Let us repeat that sequential growth means that at each tick of the clock, 
t = 1, 2, ...,T, we multiply the current N x N matrix U^^^ by an elementary s/2-block in a 
randomly chosen place i along the diagonal [see Eq. fITT]) ] with uniform distribution. The 
above analysis of "microscopic" transformation of a matrix Vn (from VN{t) to ¥^{1 + 1)) 
results in the diffusion equation for the distribution function W{fii, ...,fiN,t) 

-iy(/ii,...,/i,v,t) =/^5^ (9^7^ + '''''''"'''^') l^(/ii,...,/iN,t); (yt/7v+i = 0), (47) 

where [ij = In Aj (j = 1, A^) and Aj are the eigenvalues of a matrix V^^if). 

Suppose now that we record changes in matrix _V/v within some "mesoscopic" time in- 
terval. At, thus, producing a sequence {VAr(t = 0), yAr(At), VAr(2At), VAr(A;At)} {k > 1). 
Choosing At ~ A^, we ensure the uniform growth of a bandwidth 6 of a matrix V^r with 
time. Until the width h of the band is of order of unity, one cannot say anything concerning 
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the distribution of eigenvalues of VAr(t). Meanwhile, it is known (see |43|) that for b ^ \/iV 
one can apply the results of the Dorokhov-Mello-Pereyra-Kumar (DMPK) theory ji^] to 
the band matrices. Recall that DMPK theory describes the transmission in a system of N 
parallel interacting disordered channels. On a formal level this brings us into the subject 
of probability distribution over time-ordered products of random matrices of size N x N 
belonging to some noncommutative noncompact group G. This "dynamic" approach deals 
with the construction of a Fokker-Planck equation for the eigenvalue distribution of random 



walks on the group G (45l . |46| . The DMPK Fokker-Planck equation can be regarded as a 



heat equation for a radial part of Laplacian acting in the negative curvature symmetric coset 



space X = G/K, where K is the maximal compact subgroup of G (47|, |48|, |50 



The standard procedure implies the following steps. Introduce the radial coordinate Xi 
(system of roots) in the Cartan decomposition of X~. Note that all Xi {i = 1,...,N) are 
ordered: Xi < X2 < ■■■ < xn- The distribution of Xi {i = 1, ...A^) satisfies the heat equation 

|l^({x},t)=DA«iy({x},t), (48) 

(r) 

where D is the diffusion coefficient and Ag^ is the radial part of Beltrami-Laplace operator 

(r) 

in the symmetric space X~ . The explicit form of Ag^ acting in X~ is well known - see, for 



example, [48|, |49|, ISlI- Define by i? = {a} the root space of the Lie algebra A^. Let nia be 
the multiplicity of the root a and consider the case = (3 (for short roots) and = 1 
(for long roots). We denote by W~ the subspace of roots in R that is positive with respect 
to the Weyl chamber of R. This allows one to write Abl as 

k=l 

where 

e{x} = Y[ (sinhx^)'^'^ = Y[ |sinh^Xj - sinh^a;^!^^^ JJ | sinh2a;i|^/l (50) 

a6-R+ i<j i 

One sees that the expression fH9|) for = 1 gives the radial part of Beltrami-Laplace 
operator in the Lobachevsky plane if2, where the single root, /ii of the Cartan algebra is 
the geodesic distance on a pseudosphere. Recall that for = 1 the Liouville-type equation 
f l2^ coincides with the random walk in Lobachevsky plane. 

The substitution W^(/ii, ^n, t) = ^{fii, /iAr)\E'(/ii, fiN, t) niaps the equations (jH]), 
f l49|) onto the Schrodinger equation 



- ^^{fi,t) = H^i^,,t), (51) 



dt 



with H of hyperbolic Calogero-Sutherland-type form 



2 '^\dx^j sinh^ 2xj J 4 Vsinh^(a;fc - xj) sinh^(xA: + Xj 

where c is some irrelevant constant. 



+ —r^, r 1 + c, 

(52) 
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It has been pointed out in 5l| that the identification of fl49p - fl52p with the original DMPK 
equation |45| is possible if and only if the coupling constant in the Calogero-Sutherland 
Hamiltonian, g'^ = equals to |mQ,(mQ, — 2)|ap, where |a| is the root length. While we 

remind that nia = /3 (for short roots), rria = 1 (for long roots), in our further consideration 
we shall not fix specific value of /3, but consider it as an arbitrary parameter. This enables 
us to establish the connection between hyperbolic Calogero-Sutherland (l5T]) - (l52l) and Toda 
(1471) in the so-called Inozemtsev limit (Hi. 



To take the Inozemtsev limit in f l52l) . we should rescale the constant (3 simultaneously 
changing of coordinates Xj (j = 1, ...,N): 



/3 



(A oo) 



(53) 



Applying fl53l) to fl52|) and denoting the limiting Hamiltonian by H, we get 



2H 



Ke 



2{xj-Xj+i) 



(54) 



which coincides with the Hamiltonian of Toda chain. 

Let us end up this Section by highlight the key observation, i) On one hand, the transmis- 
sion in multi-channel disordered wires, which corresponds to the Calogero-type system for a 
special choice of coupling constant, can be described by diffusion in a symmetric space, ii) On 
the other hand, the diffusion equation in the noncompact space if at = SL{N, R)/ SO{N) 
corresponds to Toda-type system and appears in the consideration of sequential uniform 
ballistic deposition in multi-column box. The Inozemtsev limit 52| allows to establish the 
connection between these two different problems i) and ii). 



IV. TWO VIEWS ON LAX REPRESENTATION FOR PERIODIC TODA 

CHAIN 



In this Section we shall briefly consider the NNN ballistic deposition problem with peri- 
odic boundary conditions. We have already argued before that the problem is described by 
the Toda chain system. However we did not pay much attention on the boundary conditions. 
This is an important issue since the periodic and open Toda chains behave quite differently. 
It is known that the open chain is related to the motion at SL{N, R)/SO{N) while the peri- 
odic chain is related to the motion at the homogenous space of the affine algebra. Moreover 
the classical motion in the periodic case is described in terms of the Riemann surfaces whose 
Jacobians are identified with the complex Liouville tori of the dynamical system. The solu- 
tion to the classical equation in the periodic case is described via the elliptic functions at the 
Riemann surface while the quantum problem responsible for the NNN ballistic deposition 
problem is quite involved. 

Let us make a couple of general comments concerning the Toda realization of the BD 
model discussed above with the periodic boundary conditions. First, consider the classical 
Toda model which allows two different Lax representations. If there are degrees of freedom 



Random ballistic growth and diffusion in symmetric spaces 19 



corresponding to the motion at SL{N, R)/ SO{N), the Lax operator can be written either 
as 2 X 2 or iV X matrix. In the first representation the Lax operator in terms of the 
conventional phase space coordinates reads as 



^Toda 



-e-1 



(55) 



This Lax operator acts as the shift operator on the Baker-Akhiezer (BA) function along 
the sites of the chain. To compare with the discussion above consider the two-component 

Baker-Akhiezer function = ( ) . Then the linear problem can be rewritten as 



Xn 

tpn+1 - Pn^n + c''"-'^-' = Xipn, Xn = -c'^-'^l^n-l- (56) 

Putting Xn = const and freezing the dynamical variables in the system we get the equation 
equivalent to ([8]). Thus, we come to the following important conclusion: 

Statement. The first component of the Baker-Akhiezer function \E'„ plays the role of the 
finite-temperature partition function of a directed polymer in the random environment. The 
local heights can be extracted from \E^„ as in ^ by taking logarithm and passing to zero- 
temperature limit. 

This statement approves the link between the NNN ballistic deposition and the "soft" 
deposition conjectured in the Section III U[ 

Note that the 2x2 Lax operator plays the role of the Hamiltonian for the discrete 
fermionic system while (|5]) can be identified with the transfer-matrix of Toda model. 

The solution to the classical equations of motion is described by the spectral curve 

N 

det [T(A) - c^] = 0; T(A) = 1[L,{X), (57) 



i=l 



where variable u corresponds to the quasimomentum for periodic boundary conditions. 
Generically it is genus (A^ — 1) curve which gets degenerated into the rational one for the 
non-periodic case. This spectral curve is the Liouville torus for the dynamical system whose 
degrees of freedom are localized on this curve. In terms of the spectral curve one can define 
the separated variables that constitute the set of identical noninteracting degrees of freedom. 
These degrees of freedom are identified with poles of the BA function on the spectral curve. 

The question deserving special attention concerns the possibility to find a Douglas- 



Kazakov phase transition [53| in a Toda chain for different boundary condition since this 
transition is responsible for Kardar-Parisi-Zhang scaling (see js^l for details). We plan 
to discuss this question in a separate publication and as a preliminary remark we outline 
in Appendix [B] the connection between solutions of ID Toda lattice equation and rational 
solutions of Painleve II. 
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The Toda system can be also written starting from the hnear problem 

/ -pi 63(92-91) 



62 
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(58) 



with the A^-component Baker-Akhiezer function $ = {e '^"^'^^n}- This leads us to the 
spectral curve 

det {C{w) - A) = 0, (59) 



which is equivalent to the spectral curve in the 2x2 representation and reads as 



00 + UJ 



-1 



(60) 



The symmetry of two different representations of Lax operator by 2 x 2 or x matrices 
can be interpreted as "time" -like (a) and "space" -like (b) representations schematically 
depicted in Fig. [51 Geometrically this symmetry means the rotation of the picture in the 
homogeneous space - see j55| for details. In terms of ballistic deposition the time-like 
representation can be visualized as growth along the time running since the first deposition 
event, while the space-like representation can be associated with growth "along a surface" 
of the heap (i.e. in the direction orthogonal to the time axis). 



(a) 



(b) 




Figure 5. Two possibilities of constructing Lax representation for Toda chain: (a) "time-like" and 
(b) "space-like". The visualization of these representations in terms of growing heap is shown 
nearby. 



Let us conclude the reminder about two types of Lax representations by noting that 
ideologically similar description is allowed for the one-dimensional Anderson model. Two 
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different points of view on this problem, described below, deal nothing with integrability, 
however they have exactly the sarne physical interpretation as two different Lax representa- 
tions of ID Toda chain. Following [56|, consider the electron hopping on the one-dimensional 
lattice described by a tight-binding one-dimensional Schrodinger equation with site disor- 
der. Denote by ipj and Uj the wave function and the random potential on the lattice site j. 
Then the eigenvalue problem is given by the equation 



(61) 



Equation (1^ admits two representations. Introducing the two-component vector {'ipj+iy i^j)^ 
we can rewrite (161 p as 



i'j 



E 







(62) 



which has similarity with f l58|) . From the other hand, Eq. fl6T|) has representation in terms of 
spectral problem for x tridiagonal matrix C applied to the vector \& = {ipi, ip2, ■■■4'n}- 
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(63) 



It might be interesting to check whether such a structural analogy between two Lax repre- 
sentations of Toda chain and two descriptions of Anderson model has more deep physical 
grounds than just a formal similarity. 



V. CONCLUSION 



In this paper we have reformulated the NNN ballistic deposition process in terms of time- 
ordered product of the generators of a local group L^r (see the Definition 1). The generators 
of this group are the A^ x A^-matrices consisting of a single s/2^block randomly located along 
the diagonal. We have related the Markov dynamics on L^r with the diffusion in the metric 
space Hn = SL{N, R)/SO{N). 

Let us summarize the main results of our work: 

• By constructing the Beltrami-Laplace operator on Hjsf in an appropriate system of 
coordinates obtained via Gauss decomposition, we claim the existence of the integrable 
structure behind the NNN deposition problem. 

• We demonstrate the relation between the NNN growth and dynamics of the ID Toda 
chain where the number of sites in the chain coincides with the number of the sites in 
the deposition problem. The largest geodesic distance in the corresponding symmetric 
space if AT = SL{N, R)/ SO{N) is identified with the maximal local heights in the 
deposition process. 
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• Taking the Inozemtsev limit in the hyperbohc Calogero-Sutherland model we uncover 
the connection between the ballistic deposition process and the transmission in multi- 
channel random wires. This connection certainly deserves further study. 



It is worth pointing out the relation between the fluctuations of the partition function of 
a directed polymer at finite temperatures and the quantum Toda lattice has been obtained 
by N. O'Connell in a model of a semi-continuous directed polymer model with the length 



related to increments of Brownian motions [57[ • Within this model it was obtained that the 
free energy, as a function of t, behaves like the first component of a process in which 
also satisfies a diffusion equation with a Toda Hamiltonian. The extension of this work to 
statistics of vicious random walks with "a long-ranged killing term" has been considered 



in a recent work [58|. It would be very interesting to establish correspondence between the 



subject of our work, the problems analyzed in [57|, |58| and eventually to the Yang-Mills 



formulation of the statistics of vicious random walks developed recently in [5 

Let us outline the possible directions for further investigation. It seems to us that the 
most interesting case corresponds to the periodic boundary conditions. In this case the 
deposition problem gets mapped onto the periodic Toda chain whose solution to the classical 
equation of motion is described by the Riemann surface. In this case one could define the 
prepotential which yields the solution of the dispersionless Whitham hierarchy similar to 
the one considered in jHof. We plan to discuss the periodic case separately. 

Note the analogy between the ballistic deposition and the Laplacian growth (LG). In 
LG case the integrable system is defined as well, where the natural degrees of freedom are 
attributed to the shape of the droplet while the Whitham times get identified with the 
momenta of the droplet. The initial evolution equation selects the particular solution of the 
Whitham hierarchy. In the BD case the local heights play the role of the Toda degrees of 
freedom. From the physical point of view the difference between the LG and BD consists only 
in deposition rules: in Laplacian growth the particles follow Brownian trajectories until they 
are attached to the cluster, while in ballistic depositions the particles follow straight parallel 
trajectories as in a rain. The passage from LG to BD consists in "elongation" of Brownian 
trajectories in an external field. This point of view is the subject of our forthcoming work. 



Our final remark concerns the existence of shocks in Toda chain [65|. In [66[ we have 
demonstrated that the discrete stochastic equation describing the BD process in the NNN 
model can be naturally represented in terms of a dynamic programming language. This 
dynamic programming point of view allowed for a systematic translation of the study of 
growing clusters and channels separating them ( "crevisses" ) into the language of maximiz- 
ers and shocks in discrete equations of Burgers or Hamilton- Jacobi type. Whether these 
geometrical properties of a growing heap have interpretation in terms of shocks in Toda 
chain, is an open and challenging question. 
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Appendix A: The connection between the columnar heights and coordinates in the 

symmetric space 



1. We expect for > 2 the following linear relation between heights {hi, h^} and 
geodesic distances {/ii, ...,fiN} to hold: 



max{/ii, ...,h]\f} = max{/ii, /iat} 



(Al) 



where is some constant, fij = InAj (j = 1,...,N) and are the eigenvalues of the 
matrix Vt, which is the product of generators of locally free semigroup F^. The behavior 
flAip is supported by numerical simulations of NNN uniform ballistic deposition in a A^- 
columnar bounding box, considered as: i) sequential adding of T elementary blocks with 
NNN interactions, and ii) time-ordered multiplication of T matrices-generators of the locally 
free semigroup F^. 

In the Fig. [6] (a) we plot the ratio Jn(T) = max{/ii, /itv}/ max{/ii, /xat} in the 
interval T = 1000 -^ 10 000 for fixed values of (for better visibility the coordinates 
(7Ar, i/T) are used). We see that for each the coefficient 7Ar(l/T) tends for T — )> oo to 
some selfaveraged value 7Ar(0). 

The figure Fig. E] (b) demonstrates the dependence of 7Ar(0) on A^ for A^ = 2, 4, 6, 8, 
10, 20, 30, 40. The value 7Ar(0) saturates very rapidly with increasing of A^ and reaches the 
limiting value 0.8. 



0,84 




0,0000 0,0002 0,0004 0,0006 0,0008 0,0010 

1/T 



(a) 



1,0 



0,9 



-1 — ' — I — ' — r- 



0,6 



0,5' ■ ' 



4 



(b) 



_I 1 I 1 L_ 



12 16 20 24 28 32 36 40 

N 



Figure 6. Linear dependence of maximal height for NNN-deposition on the logarithm of the 
maximal eigenvalue in corresponding matrix representation: a) dependence of on 1 /T for fixed 
N; b) dependence of 7Ar(0) on different A^. 

The results of simulation depicted in Fig. E] demonstrate the linear dependence of the 
maximal height on the logarithm of the largest eigenvalue in associated matrix representa- 
tion. Moreover, for A^ ^ 1 the coefficient 7Ar becomes independent on A^. 



2. The relation between the height in A^ = 2-ballistic deposition and the corresponding 
hyperbohc distance in H2 = SL{2, R)/ S0{2) can be easily estabhshed. To see that, note 
that the height in the column is in one-to-one correspondence with the shortest distance on 
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the corresponding Cayley graph - see Fig. |l]of the Section Hi C[ Let us discuss in details the 
two-columnar deposition and show that the maximal height in two neighboring columns, i.e. 
the average number of steps along a Cayley graph is linearly proportional to the hyperbolic 
distance, d, in the space H2. We use some facts about the isometric embedding of graphs in 
hyperbolic planes and uniform Markov processes in hyperbolic domains - for more details 
see, for example, (67l |. 

Consider the hyperbolic Poincare upper half-plane H2 = {z & C, lm{z) > 0}. The 
matrix representation of the generators of the group r2 is well known (see, for example j40l|). 
however for our purposes it is more convenient to take a basis consisting of the composition 
of standard fractional-linear transform and complex conjugacy. Namely, denoting by z the 
complex conjugate of z, we consider the following action 



A possible set of generators is then: 



a b 
c d 



hi 



z 



az 



cz + d 



1 2/Vs\ , 

-1 r 







1/v^ 




(A2) 



(A3) 



Choosing the point [xq, iyo) = (0, i) as the tree root — see figJZl any vertex on the tree is 

n 

associated with an element Vn = Y\ where G {0, 1, 2} and is parameterized by its 

fc=i 

complex coordinates Zn = Ki((— l)"z) in the hyperbohc plane. 



1.4 



1.2 




Figure 7. Poincare hyperbolic upper half-plane i?2- A and B are two neighboring vertices of 
Cayley tree isometrically embedded into H2; d is the hyperbolic distance between points A and B. 

For H2 = SL{2, R)/ S0{2) one can identify an element with its class of equivalence of 
S0{2). If one denotes by fi{Vn) = fi{i,Zn) the hyperbolic distance between i and z^, the 
following identity holds 

2cosh/.(K) = Tr(KKl) (A4) 
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where dagger denotes transposition. 

Let WnifJ') be the probabihty to find the tree vertex in the generation n at the hyperbohc 
distance fi from the root point. In terms of deposition process, n is the average height of 
two neighboring columns). It means that we are looking for the distribution of the traces for 
matrices Vn which are irreducible products of n sequential multiplications. If we denote by 
the irreducible length of the word represented by the matrix Vn, then Vn is irreducible 
if and only if fi{Vn) = n. Such a word enumeration problem is simple in case of the group 

n 

T2, because of its free product structure: T Z2 ® Z2 ® Z2. Indeed, if Vn = JJ^ ha^ one 

k=l 

has iiiVn) = n ii and only if ha^ 7^ V/c. Hence we have to study the behavior of the 

random matrix Vn-, generated by the following Markovian process 

+ 1) mods with probability \ 

(A5) 

(a„ + 2) mod 3 with probability | 

One can use the standard methods of random matrices and consider the entries of the 
2 X 2-matrix Vn as a, 4-vector Vn- The transformation Vn+i = Vnha reads 

K« = ( 0' "t ) V. (A6) 

This block-diagonal form allows one to restrict ourselves to the study of one of two 2- 
vectors, composing V„, say Vn- Parameterizing Vn = {gn cos 9n, Qn sin 9n) and using the 
relation = /i„ ^ 21n^„ valid for n ^ 1, one gets a recursion relation Vn+i = hl^Vn in 

terms of hyperbolic distance fin- 



fin+l = fJ-n + In 



5 4 

- + - cos(26'„ + (fc 



(A7) 



where ipa depends on the transform through ipa = {2a — l)7r/3 (a = 0, 1,2), while for 
the angles one gets straightforwardly 

tan6'„+i = ha(tan9n) (A8) 

and the change a ^ a means the substitution (0,1,2) — (1,0,2). Action of ha is still 
fractional-linear. The properties of the invariant measure z/q(^) have been discussed by 



Gutzwiller and Mandelbrot [68|. We are now led to compute the Lyapunov exponent, 7 
which is 

^ = lim M = f' dO^iie - ^) In f I + ^ cos 2^^ (A9) 

n-!>oo n Jq b \6 6 J 

The numerical value of 7 can be obtained by means of semi-numerical pro cedure which 
involves the numerical information about the invariant measure /ii(^) [g^, 6^- One finally 



gets: 7 ~ 0.79. By this we prove the linear relation = 777, between hyperbolic distance, n, 
in H2 and average height in two neighboring columns, n. 
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Appendix B: Minimal action principle for Toda system and rational solutions of 

Painleve II 



Represent the partition function W in the form of Feynmann-Kac path integral with the 
Lagrangian obtained by Legendre transformation from Toda Hamiltonian 



"Dj/i} exp 



ds 



■ N 

^fi^s) + v(^fij{s),fij+i{s] 



, (Bl) 



where V^(/^j, /Xj+i) = e'^-' '^^+1 is the Toda potential Minimizing the corresponding action, 
we arrive at the standard Euler equations, which describe the dynamics of ID Toda chain: 



(B2) 



where jlj{s) = -^jjj{s). There are many equivalent representations of Toda chain equation 
(lB2p . Among the most known, is the bilinear form of (lB2p (see, for example, 60|): 



'j 'J 



where = f r,-, r/ = and 



fij = In 



(B3) 
(B4) 



The tau-function r„ can be written in the determinant form (see, for example, [6l| 



Tj = det 



d \2i-2 



(B5) 



The function ip for open boundary conditions is an arbitrary function of s. Applying a gauge 
transformation to To 



with arbitrary p and g, one arrives at the equation 



2 In Tj (s) = In cTj (s) + pg. 



( l\2 2 



3 ' 



(B6) 



(B7) 



where (T_i = ctq = 1, (Ji = g. Equation flB7p resembles the differential-difference equation 
that defines Yablonskii-Vorob'ev polynomials j62[, Qn{z): 



Q''{z)Q,{z)-{Q'{z)f 



4 



(B8) 



with Qo = 1, Qi = z and the derivatives taken with respect to z. To establish the precise 
connection between (IB7p and fIBSp . let us make the shift in f lB7p : aj — )■ A-' aj and then fix 
the values of A, p, and q to satisfy conditions: A = \, pq = — |. 
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It is known 63|, |6J] that the Yablonskii-Vorob'ev polynomials Qni^) are connected to the 
rational solutions of Painleve II. Namely, the combination 



^ dz Qj{z) - Q,_i(^) Q,{z) 
gives the rational solutions of the Painleve II equation 

w"{z) = 2w^{z) + zw{z) + a {a = jeZ). 



(B9) 



(BIO) 



So, the classical optimal trajectory solution for Toda Hamiltonian can be linked by equations 
(IBip - flBlOp to the rational solutions of the Painleve II equation. 
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